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We generalise the relativistic expression of Ohm's law by studying a multi-fluid system of charged 
species using the 1+3 covariant formulation of general relativistic electrodynamics. This is done 
by providing a fully relativistic, fully nonlinear propagation equation for the spatial component 
of the electric 4-current. Our analysis proceeds along the lines of the non-relativistic studies and 
extends previous relativistic work on cold plasmas. Exploiting the compactness and transparency 
of the covariant formalism, we provide a direct comparison with the standard Newtonian versions 
of Ohm's law and identify the relativistic corrections in an unambiguous way. The generalised 
expression of Ohm's law is initially given relative to an arbitrary observer and for a multi-component 
relativistic charged medium. Then, the law is written with respect to the Eckart frame and for a 
hot two-fluid plasma with zero total charge. Finally, we apply our analysis to a cold proton-electron 
plasma and recover the well known magnetohydrodynamic expressions. In every step, we discuss the 
approximations made and identify familiar effects, like the Biermann-battery and the Ifall effects. 

I. INTRODUCTION 

Relativistic plasmas play a major role in high-energy phenomena, such as those associated with active galactic 
nuclei, black-hole magnetospheres, relativistic jets, the early universe, etc. Although a complete kinetic description 
of a relativistic plasma is desirable, it becomes less efficient when the large-scale bulk plasma motions are dominant. 
Thus, multi-fluid hydrodynamics is more appropriate in astrophysical systems characterized by those kind of motions, 
and in this framework it is desirable to coarsen the description to that of a one-fluid magnetohydrodynamics (MHD). 
The MHD equations, which comprise of the particle conservation law, the energy/momentum conservation laws and 
Maxwell's formulae, must be also complemented by Ohm's law, namely an equation relating the induced electric 
current with the electric field of the plasma. Ohm's law is directly involved in the magnetic induction equation, used 
in the description of bulk plasma dynamics. Despite its importance, however, there is still no general consensus on 
the form of Ohm's law for relativistic plasmas. 

Some of the previous work towards obtaining a relativistic version of the generalised Ohm's law can be found 
in [U, 0, H, 0, H, [1] . In general, the authors start from a kinetic description and subsequently derive the required law 
in the form of an electric-current conservation equation. The approach of ^ assumes a multifluid description, but the 
study is then confined to a specific black-hole metric. Here, we will derive the relativistic version of Ohm's law for a 
multi-component plasma, a general metric and in the presence of an arbitrary electromagnetic field. We will do so by 
means of the 1-1-3 covariant approach to General Relativity and Electromagnetism (see 0, H, S [l0| for reviews and 
further references [2l|). Ohm's law is given in the form of an evolution equation, along a Fermi-propagated frame, for 
the spacelike component of the electric 4-current. Exploiting the mathematical compactness and physical transparency 
of the covariant formalism, we arrive at an expression that allows for a direct comparison with the existing Newtonian 
version of Ohm's law and identifies the relativistic corrections in an unambiguous way. An additional advantage is 
that the covariant version of Ohm's law provided here goes beyond the cold plasma limit of the expressions obtained 
so far and applies to fully relativistic fluids. 

The aim of the present, first, paper is primarily to set up the mathematical framework, discuss the physics and 
provide the relativistic expression of Ohm's law in its most general possible form. The latter means that our results can 
be applied to a wide variety of situations, ranging from relativistic plasmas to astrophysics, large-scale structure and 
cosmology. Thus, given any arbitrary "background" metric (astrophysical or cosmological), one can easily linearise our 
the equations around it. In addition, the H-3 covariant formulation is built on irreducible kinematical and dynamical 
quantities, which assigns a clear physical interpretation to every variable in our equations. The latter should prove 
very advantageous when studying nonlinear effects, particularly those outside the ideal-MHD limit, like current-sheet 
formation, turbulent plasmas, magnetic-dynamo amplification and dissipative effects. Applications of this kind will 
be the subject of future work. Here, maintaining the generality of the metric and of the electromagnetic field, we 
focus on a two-component fluid and apply our multi-fluid expression of Ohm's law to a system of relativistic protons 
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and electrons. Adopting the Eckart frame (e.g. see [l^), allows us to define the bulk quantities, such as the plasma 
velocity for example, in a way that closely resembles those of the non-relativistic treatments. This, in turn, helps us to 
discuss the common approximations in the use of the law and to explicitly show the relativistic counterparts of very 
well known Newtonian effects, like the Biermann battery and the Hall effect. 



II. MULTI-COMPONENT FLUIDS 



A. 1+3 kinematics 



Consider a pseudo-Riemannian spacetime filled with metric gat of signature (+,+,+,—) and introduce a timelike 
velocity field Ua normalized so that Wqu" = — 1. The tensor 

hab = gab + UaUb , (1) 

projects orthogonal to Ua and into the observers local rest-space. [12] The Ma-field defines the frame of our fundamental 
observers and, together with hat, achieves a unique 1-1-3 "threading" of the spacetime into time and space. The same 
two fields are also used to define the covariant time ans spatial derivatives of any tensor field Sab - '^'^ according to 

Sab-''' - = u'VeSab-''' - and B.Sab-'^ - = h/ hj hg'h/ S fp...'^'- ' ' , (2) 

respectively (see @, [13] for reviews on the covariant formalism and more details). 

Covariantly, the kinematics of the fundamental observers are determined by decomposing the orthogonally projected 
gradient of their 4-velocity field into its irreducible parts. In particular, we have 

DfcUa = ^ Qhab + (Tab + t^ab " UaUb , (3) 

where Q = V°Ua = D^Ua determines the average expansion (when positive) or contraction (when negative) of the 
volume element associated with Ua, (Jab = 0(6^^^ indicates changes in its shape (under constant volume), LOab = D[(,Ua] 
measures the rotational behaviour of Ua and iia — u^^bUa is the 4-acceleration. The latter indicates the presence of 
non-gravitational forces. 



B. Tilted 4- velocities 

Let us assume a mixture of fiuids and define the 4-velocity of the i-th fiuid component (with u'^^^Ua ^ — —1) as 
the future directed timelike eigenvector of the associated Ricci tensor. Then, the tensor 

/i^i^=. 9-^^ + 4)4) (4) 

projects orthogonal to The relation between and u"^^ is determined by the hyperbolic angle /3(j) between the 
two 4-velocity vectors. This 'tilt' angle, which also determines the peculiar motion of the i-th fiuid component relative 
to (see Eq. © below), is given by [13l] 

cosh/3(j) = -UaU^i-^ , (5) 

with 

4)=7(o(^'' + ^w) (6) 

and w^j-j representing the peculiar velocity of the i-th species relative to the Ua-frame. Also, 7(i) = (1 — f^^^)"^/^ is 

the Lorentz-boost factor and v'^-^ = v'^--jVa^ . When the tilt angle is small (i.e. for f3^i) <ti 1) we have ~ This 
means non-relativistic peculiar velocities and 7(i) ~ 1. 
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C. Multi-component perfect fluids 



A general material medium is described by its energy momentum tensor T°^, its particle flux iV" and by the entropy 
flux vector S"'. The former two quantities are conserved (i.e., VbT"''^ = = Va-/V°), while the last obeys the second 
law of thermodynamics (i.e., Va<5'° > 0). When the strong energy condition holds, the energy momentum tensor of 
a fluid accepts a unique timelike eigenvector u?p, normalised so that u^u^ = — 1 UM- One may also define a unitary 
timelike vector parallel to N°- by = N'^/ ^/—NaN"-. Provided that the fluid is perfect ( or in equilibrium), f/^ ^, u^ t 
and S*" are parallel and define a unique hydrodynamic 4- velocity vector, the rest-frame of the fiuid flow (see e.g. |17|'). 
This is also the only frame in which the energy momentum tensor of the matter assumes the perfect-fluid form. When 
dealing with an imperfect fiuid, however, there is no longer a uniquely defined hydrodynamic 4-velocity. 

Consider a mixture of perfect fluids, where the i-th component has energy density isotropic pressure and 
moves along the timelike 4-velocity field u^^-j. Relative to this frame, the energy momentum tensor of the individual 
species decomposes as 

with h'^-^ given by (|4]). Also, the associated particle fiux vector is given by 

= , (8) 

where n(i) is the number density of each matter component in their own frame. 

Substituting the transformation law ^ into Eq. ^ allows us to re-express the latter with respect to the Ua-frame. 
The result reads 

Tab ^ a I ^ l^dh i a ^fc i ^> i ^ah frk\ 

(i) = u + P{i)h -Hu +7r(i) , (9) 

with ft,"'' given in ([T]) . The above corresponds to the energy-momentum tensor of an imperfect fluid with 

Aw = {^^{^) + Pit)) ~ Piz) , (10) 

m = P(i) : (11) 

9w = 7(')(A*W+m)«W' (12) 



- 7f.)(A^(.)+Pw)«w^'^,o- (13) 



We emphasise that the forms of both and tt^^^ are different from those usually found in the standard literature 
(e.g. see [13, IS]). In particular, 7t°^^ is no longer trace-free, which means that its does not only contain the purely 
anisotropic part of the system's pressure. For the same reason pj-j) is not the effective isotropic pressure of the total 
medium, since part of the latter is incorporated in the 7rj?^^-tensor. Neverthelees, the adopted representation has serious 
technical advantages, which will allow us to go easily beyond the cold-plasma limit without compromising the physics 
of our discussion. 

As with the energy-momentum tensor, we may also express the particle flux of the species relative to the it°-frame. 
To be precise, combining Eqs. ^ and ^ we arrive at 

N^.-^^h^^u'^+U^^, (14) 

with 

"W=7W'^W and Sf^^ = fn^^v"^-^ , (15) 

representing the associated number density and particle flux respectively. When dealing with non-relativistic species, 
7(j') ~ 1 and we may ignore terms of quadratic (and higher) order in u^^j. Then, according to (|10p - (ll3p . we have 

A(i) - /^(i): P(i) - P{i), - (/^(i) +P{i))'"{i) and tt^'] ~ 0. Also, for vf-^ < 1, relations ^ reduce to h(i) ~ and 
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III. MULTI-COMPONENT CHARGED FLUIDS 



A. Electromagnetic fields 

The Maxwell field is covariantly characterised by the antisymmetric electromagnetic (Faraday) tensor F""^. Relative 
to an observer moving with 4-velocity u", the later decomposes as (e.g. see 0,0)) 

F""^ = 2u[°£;''l + e"''"Sc , (16) 

where e"^"^ is the permutation tensor orthogonal to u° and E"- — F°'^ui,, = e'^''^Ff,c/2 are respectively the electric 
and magnetic field measured by the fiducial observer. The evolution of these two fields is monitored by means of 
Maxwell's equations, which in covariant form read 0, @ 

^(a) ^ _ 2 ^ ^^^^ ^ uj''h)E^ + e^bcU^B" + curlS" - J" , (17) 
3 

B^""^ = 95° + (cr'^b + t^%)B^ - e'^bcU^E'^ - cuiW , (18) 

DaE" =p- 2uJaB'' (19) 

and 

D,S" = 2L0aE'' . (20) 

Note that J'^ = h°-bJ^ and p — —UaJ"" are the spatial current and the charge density respectively, with Ja representing 
the 4-current.[23l| Also, uj"' ~ e°'bc^^'^ 1^ is the vorticity vector (with — e"''cW'^) and curlw" = e°'bJ^^v'^ for any 
orthogonally projected vector v°' . 

The Faraday tensor obeys Maxwell's equations and also determines the energy-momentum tensor of the electromag- 
netic field by means of the familiar formula 

T^L) = F^'F'c - \ F-^F^ag'^'' . (21) 

(em) 

Maxwell field. Thus, relative to the Ua-frame, 

T(L) ^\{F'^ + B^) + + B^) + '^.Q^'^u^^ + , (22) 

with £'2 = E^Ea, S2 ^ B'^Ba, Q" = e'^^'^EbBc and = -E^'^E''^ - B^"-B''K In other words, the electromagnetic 
field can be treated as an imperfect medium with energy density (E'^ -|- B^) /2, isotropic pressure (i?^ + i?^) /2, an 
energy flux represented by the Poynting vector and anisotropic stresses given by 11°'' 0, Q . 



The latter combines with Eq. (fT6|) to facilitate the irreducible decomposition of TT!.^ and a fluid description of the 



B. Conservation laws 



Consider a multi-component system containing species of different nature (e.g. baryonic and non-baryonic matter, 
photons, etc) in the presence of an electromagnetic field. The charged particles are coupled to the Maxwell field and 
the mixture has a total energy- momentum tensor given by the sum T^^ + T^^^^ . The latter satisfies a conservation 
law of the form 

V,T(°^ - F\Jl) = , (23) 

since '^bT^^^) ~ ~F"'bJ^iy where J^"^ is the electric 4-current density of the i-th (charged) species. The interaction 
term Q^^-^ represents forces other than electromagnetic. The latter, as a consequence of the conservation of the total 
energy-momentum tensor, obey the constraint 



i 



(24) 
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The timelike and spacelike parts of expression (j23p respectively provide the conservation laws of the energy density 
and of the momentum density of the i-th fluid component. Thus, by projecting ((23)) along we arrive at the covariant 
form of the generalised continuity equation 

ki) = - + m) - Dag",) - 2Uaq1^) - (^<Jab + ^ Qhab^ Trg + EaJ^^-, ~ G^^) , (25) 

where J7^°) = h°'bJ^i) and Qi^i) = UaG^^iy Note the second last term in the right-hand side of the above, which describes 
alternations in the energy density of the i-th fluid due to the action of the electromagnetic field. This term may be seen 
as representing the familiar Joule heating effect in a covariant manner. Recall that tt"^^ is not generally a trace-free 

tensor and therefore the sum habT^^i-^ does not necessarily vanish. Similarly, projecting Eq. orthogonal to it", gives 
the covariant form of the generalised Navier-Stokes equation 

+ (p(,)£;'^ + e%e:^('^)i3^)+4f , (26) 

with q"^}^ ~ ^°b9(i) and = —UaJ^^-^ representing the charge density of the i-th component. 

Turning to the particle numbers and assuming that the individual number densities are not necessarily conserved, 
we may write 

^aN^^ = Q(,) , (27) 

where Q(i) indicates either an increase or a reduction. Particle-antiparticle annihilation, for example, will reduce the 
numbers, while = ensures overall particle conservation. Substituting (fT4|) into the above and then using 

expressions leads to 

n{i) = - {q + DaUf'j)) - W('j)DQn(i) - UaW('j)n(i) -f , (28) 

thus providing the propagation equation of the particle's number density relative to Ua- We close this section by 
pointing out that both Eq. ([^5]) and Eq. (pS)) contain implicit derivatives (temporal and spatial) of 7(,;) . This is the 
consequence of our decomposition choice (reflected in the set ([^- (|13p ) and ensures that we can treat relativistic (hot) 
plasmas with minimal technical complexity. 



C. Electric charges and currents 

The 4-current density of each charged species is related to its associated particle- flux vector by J"^.-^ — eZ(^i-jN^-y 
with e representing the fundamental electric charge and the atomic number of the particles. Recalling that 
-^{i) ~ ^^^H'^'^ + '^{i))^ relative to the m° frame - see Eqs. (ITil) . ([TS|) . we have 

J^,^^eZf^,)h^,^{u^ + v^^^^) . (29) 

The timelike part of the above gives the charge density of the corresponding charged component in the fundamental 
frame, while its spacelike counterpart leads to the associated 3-current. In particular, projecting (|29p along and 
orthogonal to Ua we arrive at 

P(i) = -UaJ^,) = eZ(,)n(,) and J('J) ^ j'l^^ = eZ(,)n(,)W^,) , (30) 

respectively. Consequently, the total charge and the total 3-current are given by the sums 

P^^Pi^)=^^^i^)Hi) and '^'^ = Xl'^W = ' (^1) 

i i i i 

with the former vanishing in the case of overall electrical neutrality. The latter is a good approximation on scales larger 
than the Debye length of the species, where the bulk properties of the plasma dominate. In that case = 0. 
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IV. GENERALISED RELATIVISTIC OHM'S LAW 



Relativistic expressions of the generalised Ohm's law, in the form of a propagation equation for the electric 3-current, 
appear in various versions in the literature [JS, 0, With the exception of however, all of the aforementioned 
studies address cold plasmas only. Here, we use the 1+3 covariant approach to relativistic hydrodynamics and electro- 
dynamics to express the generalised Ohm's law in terms of irreducible kinematical, dynamical and electrodynamical 
variables. The result is an expression that closely resembles the non-relativistic forms encountered in the standard 
plasma literature on one hand, while on the other it identifies the relativistic corrections in an unambiguous way. 

Covariantly, the time evolution of any given quantity is monitored by the orthogonally projected time-derivative 
of the associated variable. When the latter is spacelike, as it happens in the case of the 3-current, the orthogonally 
projected time-derivative coincides with the familiar Fermi derivative. Thus, according to (j31[) . 



J 



a) 



(32) 



To obtain the full expression, one needs to replace n(i) and u^^^. The evolution of h(i) is readily given by Eq. 
The one for wj^^^, on the other hand, is obtained by replacing (pS)) into Eq. (pS)) . Then, 



.(a) 



where M(,) = +P(i)- Substituting the above result into expression ([5^ . using ([^5]) and employing some lengthy 
but relatively straightforward algebra, we arrive at 



Hi) 



M, 



M, 




M, 



EbJ(i) 



(34) 



This monitors the evolution of the total 3-current, associated with a multi-component system of non-comoving charged 
perfect fluids. Given that Eq. (|34|) relates the total 3-current to the electric field, it also provides the fully relativistic 
version of the generalised Ohm's law. Recall that the "hat-variables" contain the full Lorentz-boost factors (i.e. the 
7(i)'s - see definitions pII)) - (IT^ and also (|15b .)). which guarantees that our analysis applies to hot as well as to cold 
plasmas. We also note that, in principle, one can invert Eq. ([M)) to an expression for the electric field. Combined with 
(fTS)) . the latter can be used to study the magnetic component of the Maxwell field in detail. 

Comparing the right-hand side of expression ([34]) to that of its Newtonian counterpart (see Eq. (3.5.9) in [23] for 
example), we can immediately see the analogies and also locate the relativistic corrections. Thus, the first two terms in 
the right-hand side of the above are due to the relative motion of the fundamental observers. Both terms also appear 
in Newtonian treatments and represent changes in the 3-current density triggered by inertial forces. The fourth term 
will lead to the familiar Hall effect (see ? IV CI below), while the fifth comes from variations in the effective 'isotropic' 
pressure of the individual species. Note that both the spatial and the temporal pressure gradients are involved, with 
the latter treated as the relativistic correction to the Newtonian Biermann battery effect. The sixth term in the right- 
hand side of (|34p vanishes when global electrical neutrality is imposed and the seventh is due to particle collisions 
that lead to momentum transfer between the species. The latter provide a measure of the electrical resistivity of the 
total medium. The first component of the eighth term comes from changes in the number density of the species due 
to particle-creation (or annihilation) processes, while the second is triggered by energy-density exchanges between the 
individual fluids and the Joule- heating effect (see Eq. ((25|l in IIII B|) . The second last term in the right-hand side of 
P4p accounts for spatial inhomogeneities in the velocities and the number densities of the individual species, while the 
last is triggered by relative motion (inertial) effects. The former has a Newtonian counterpart (e.g. see (3.5.9) in (20j). 
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V. HOT TWO-FLUID PLASMAS 



A. Eckart frame vs Landau frame 



A mixture of relativistic ideal fluids does not behave as an ideal medium, which means that when considering the 
bulk motions of the total system we must first specify the reference frame we are working in. Traditionally, two choices 
are in order: the Eckart (or particle) frame [12l | and the Landau (or energy) frame |14l |. For our multi-component 
fluid, the 4-velocity associated with the total particle flux is 

Wh = ^ = y^TT ^ +-^^(1)) - (35) 

with the second equality resulting from decomposition (|14p . On the other hand, using ([9l)- (|13p . we may write the 
4-velocity associated with the total energy flux as 



E ^W-^ = ^ E + 9w) ■ (36) 



Observers following the Eckart frame see no particle flux, while for those in Landau coordinates the energy flux 
vanishes. Thus, in our case, the Eckart and Landau frames are deflned by demanding that 

EA'w^O and E^W^O' (3^) 

i i 

respectively. The immediate consequence is that an observer in the Eckart frame detects a non-zero heat flux, while 
its Landau counterpart sees a particle drift. On these grounds, the bulk velocities of the plasma in the Eckart and the 
Landau frames are respectively defined by 

^E^^E^h -2^^E-^W' (38) 

i i 

where M = ^(i) ^"^^ = Si '^{i)- Both frames are physically equivalent and choosing one against the other is a 
decision dictated by the particulars of the problem in hand. In Appendix A we discuss in more detail the definition of 
these frames and their non-relativist limits). Here, we shall work in the Eckart frame because there the definitions of 
the bulk quantities (e.g. the bulk velocity of the plasma and its energy density) closely resemble their non-relativistic 
associates. This will simplify our analysis and facilitate its physical interpretation. 



B. Globally neutral plasmas in the Eckart frame 



So far our analysis has been completely general, covering multi-component systems with species of different nature 
(i.e. baryons, non-baryons, photons, etc). Most physical plasmas, however, are treated as two-fiuid mixtures of 
oppositely charged particles. In addition, Eqs. ([55]) can now be "inverted" to express the velocities of the individual 
species in terms of bulk variables only in two-component systems. For these reasons, and also for mathematical 
simplicity, we will from now on confine to plasmas containing one positively and one negatively charged species. All 
other possible constituents (as, e.g., photons) will be considered as external. Therefore Eq. (|24p must be rewritten as 

where g^^t represents interactions that now are external to our system as, for example, collisions with photons (i.e., 
Compton effect), anomalous resistivity due to scattering on turbulent flows, etc (see Q for a short account of the 
possible external effects). 

Assuming that and w^i are the associated velocities relative to the fundamental u'^-frame, the corresponding 
energy-flux vectors are 



(40) 
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with M± ~ fi± + p± — 7±(Ai± + P±) (see Eqs. pUjl - p^ *). Combining definition (pSb ) and expression pOj) we deduce 
that, with respect to the Eckart frame, the bulk velocity of the system is 



v'k^^ [M+vl + M^vt j , (41) 

where M — M±. Also, following Eq. (pib ). the total 3-current associated with our two-component system reads 

= e {Z+n+vl + Z^h^vt) , (42) 

with Z± representing the atomic numbers of the positive and the negative charges respectively. 

Proceeding in line with the non-relativistic studies, we will now express the velocities of the two charged components 
in terms of the bulk properties of the plasma. The first step towards this direction is to assume overall charge neutrality. 
This applies to scales larger than the Debye length of the species, where the bulk properties of the plasma dominate. 
In that case, Zj^hj^ = —Z_h_ = fiE - see definition ([?Tk ) - and expressions recast as 

p± = ±ehE , J| = ieniswi and J" = efiE « - v°_) , (43) 

respectively. Then, using (|43b) we can invert Eq. (|4H) and arrive at 

vl ^ v'^E ± J"" , (44) 

efisM 

which substituted into Eq. (143b ) gives 

A'U 

J2 = ±en£f B + ^ J" . (45) 

M+ + M_ 

These last two results express the velocities and the 3-currents of the individual species in terms of the corresponding 
bulk variables, provided global electric neutrality holds. 

Returning to the Eckart frame, we recall that — there (see definition (|37p ). This immediately implies that 

n+v°j_ = —n^vt- The latter combines with relations (f44|) and, given that the total charge density is zero, leads to 

enEiiyZj^ — Z_j 

Therefore, the bulk velocity of the plasma, as measured in the Eckart frame, is colinear to total current. Finally, 
substituting the above into Eqs. (|44l) and (l45l) leads to 

vl= . ^ . Ja. and J'±-± ^ Ja, (47) 

enE{Z+ — Z_) Z^ — z_ 

respectively. 



C. Ohm's law for hot two-fluid plasmas 



The last two relations monitor the kinematics of the species in terms of the bulk quantities and with respect to 
the Eckart (particle) frame, provided overall electrical neutrality holds. That aside, no other restriction has been 
imposed and the individual charged species are completely general and fully relativistic. In what follows we will use 
the results of § IV Bl to obtain Ohm's law for a relativistic two-fluid plasma. Thus, recalling that zero total charge 
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implies = 

j{a) 



-Z^fi^ and employing some straightforward algebra, Eq. 



recasts to 



4 



„,6 T3C 



M_ 



eriE 



M+M 



y^Tf 



Qextv'\. 



M4 



+671^ 



(48) 



For the physical interpretation of all the right-hand side terms and a comparison with their Newtonian counterparts, 
we refer the reader to the discussion given in § IIVI after Eq. (IM)) . Relative to that expression, we have kept the 
velocities of the individual species and replaced the associated 3-currents by means of (j43b ). The only exception was 
when dealing with the fourth term in the right-hand of (|34p . where the 3-currents of the species were replaced using 
relation (|45l) . The latter has allowed us to include the magnetic convection term and the Hall term in the right-hand 
side of Eq. explicitly. Keeping the fluid velocities, on the other hand, will help obtain the non-relativistic limit of 
the above (see § IVII below) . Also note that the condition of global charge neutrality has removed the 4-acceleration 

term from . Finally, for the interactions we wrote as q''^'' -\- G-'^ = G^xl , G+ + G- = Gext and assumed that 
Q+ + Q_ = 0. The first two conditions respectively guarantee momentum and energy-density conservation in absence 
of external interactions, while the last ensures that the overall particle number is preserved. 

Expression (|48p can be given in a variety of forms depending on the problem in hand. For example, using relations 
(|T7k ) and (147b ) we can recast the right-hand side of (IM)) in terms of the total 3-current. Alternatively, one can 
employ Eq. (I46|) to express everything in terms of the bulk velocity. Here we will do the former, while maintaining 
the condition of global charge neutrality. To compactify the results, we also introduce the auxiliary bulk variables 



M = M+ 

which immediately imply that 4M+Af_ 
4 



M2 



J 



Ha) 



E 



and A = Af+ - M_ 

A^. Employing these relations, Eqs. 

2enE[M{Z+ + ZJ) - A(Z. 



(49) 



transform into 



Z-)] 



2[M{p+Z^ 



(M2 - A2) 
;5_Z_)- A(j5+Z++|5_Z_)] 



Sens ^(a) 



M ^ 



(M2 

Qext 



-A2)(Z+ 

'iefiEM 



Z\ 



ja 



AP - A2 



£ bcJ O 

M (D''j5+ - D''p_) - A (D°p+ 



{NP - /S?){Z+ - Z 
2efiE 



A " M2 

Z-)Q- 



G': 



a) 



'^Z-f-Gext 



TIE 

z'i - zt 



- A2 

2[A(Z2 



[Z, 

ZD 



-Z-) 

M{Zl 



M + A 

+ zl)] 



{AP - A^){Z+ - Z^ 



Ja 



E^'JbJ" 
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M(Z+ + Z_)- 


A{Z+ - ZD 


G- 


iz+-zD 


'aP - A? 





Ja 



z. 



e(Z+ - Z_) 



■D" 



riE 



■JbJ'' 



e^hUZ+ 



Z) 



(j'^JbJc J" 



(50) 



with M2 — A2, Z+ — Z_ 7^ 0. Relations and ([50]) provide the relativistic (1-1-3 covariant) version of the generalised 
Ohm's law, with respect to the Eckart frame, when applied to a mixture of two hot and interacting charged fluids. 
We also remind the reader that these results have been obtained under the assumption of global electrical neutrality. 



D. Relativistic particle-antiparticle plasmas 

Expression (jSOp is particularly useful when dealing with particle-antiparticle pairs. An electron-positron plasma, for 
a example, has Z± = ±1 and M+ = Af_. The latter means that M ^ 2M±, A = and AP - A^ = AP = 4M^ (see 
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definitions In such an environment Eq. ([50]) simplifies to 

3 ' M M M 

- i E^'JbJ'' + (\qJ^ + ^"^JbJ^ J"" , (51) 

and contains no Hall effect (since A = in this case). Also, given that p+ = p_ for particle-antiparticle pairs, only 
part of Biermann-battery effect survives (the relativistic - carried by the fourth term in the right-hand side of the 
above). This and the Joule-heating term are the only purely relativistic corrections. 

VI. OHM'S LAW FOR COLD PROTON-ELECTRON PLASMAS 
A. Non-relativistic limit 

Cold plasmas have components with non-relativistic relative velocities. Thus, at the low velocity limit, one can ignore 
terms of quadratic (and higher) order in v'^^y As a result, 7(j) ~ 1, n(i) ~ n(i), ~ and ~ ^(i) ~^P(i) — l^ii)- 
If, in addition, the plasma is a mixture of protons and electrons, we have Z± — ±1. When the condition of overall 
electrical neutrality is also imposed, we may set n± = n and p± = ±en. Then, assuming that and to_ are the 
proton and the electron masses respectively (with ^± = nm± and to_ <C to+), we may write 

M = M+ + M_ ~ nm+ and A = M+ - Af„ ~ nm+ . (52) 

Note, however, that 

_ A^ = 4M+M_ ~ in^m+m^ . (53) 

Applying the non-relativistic limit to (I48p immediately removes the last four terms from the right-hand side of 
the latter (all quadratic in v^.). For cold species we may also ignore the p - terms, though the spatial variations 
of the pressure arc not necessarily negligible. When dealing with proton electron systems, the eighth (Q+) term in 
the right-hand side of Eq. (|48|) vanishes (this is clearly shown in ([50|) - recall that Z± = ±1). We may also ignore 
energy-density changes due to collisions between the cold components, and therefore remove the ninth term from the 
right-hand side of Finally, when applied to our two-component fluid, the coUisional (^1°^) term in which 

triggers changes in the momentum of the species, takes the more familiar form 

- + gi") . -^j^ . _,£!!1^„ , (54) 

M+M_ m_ 

where v is the ave rag e collision frequency and rj — vra- je^n is the (scalar) electrical resistivity of the two-component 
medium (e.g. see (20[). Note that, although here we have adopted the common approximation of a scalar electrical 
resistivity, our analysis also applies to general fluids with anisotropic (tensor) resistivity. This can be done by using 
kinetic theory to specify the interaction terms. On these grounds and by using the auxiliary relations (|52p . (j53p . 
expression reduces to 

6 m_ TO- TO- 

+ — D>_-r;— J^ (55) 

TO- TO- 

in agreement with the expressions found in the standard literature (e.g. compare to Eq. (3.5.9) in [20l|). Alternatively, 
we may recast the above into the more familiar form 



— e%cys= + — D>_, (56) 
en en 



which immediately shows the terms responsible for the Hall and the Biermann-battery effects - see the last two terms 
in the right-hand side. Either of these two expressions provides the l4-3 covariant form of the generalised Ohm's law 
for a cold proton-electron plasma. 
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B. Magnetohydrodynamic limits 

Adopting the commonly used magnetohydrodynamic (MHD) approximations, namely ignoring all the right-hand 
side terms in Eq. ()56|) . leads to the usual form of Ohm's law for an electrically resistive medium 

+ e\,v%B- = f]J'' . (57) 

We note here that, instead of dropping the 3-current terms seen inside the square brackets of (j56p . one could in 
principle incorporate them within the resistive term in the left-hand side of that expression. Then, we will be referring 
to a "generalised resistivity" that accounts for changes in the electrical properties of the medium due to relative motion 
(i.e. inertial) effects. It is also clear that in the idealised case of a perfectly conducting medium, namely for — > 0, 
Eq. ([57]) reduces to the ideal-MHD form of Ohm's law. The latter is given by the simple and well known formula 

E"" = -e\,v%B^ . (58) 
VII. SUMMARY 

Hot plasmas are of major importance in a variety of physical phenomena, ranging from laboratory physics to 
astrophysics and cosmology. A key factor in determining the behaviour of plasmas is their electrical properties and 
these are theoretically monitored by means of Ohm's law. The latter appears in a number of different versions, which 
depend on the specifics of the problem in hand. Here, we are providing fully relativistic and fully nonlinear expressions 
for the generalised Ohm's law for plasmas, by deriving the 1-1-3 covariant propagation equation of the 3-current density 
associated with a multi-component fluid. Adopting a suitable definition for the irreducible variables of the matter 
fields (see Eqs. (ITU1) - P^ and (fT5b )). we were able to address hot plasmas with a fluid-based approach and without the 
need of kinetic theory. The use of the covariant methods has also facilitated a mathematically compact and physically 
transparent presentation of the subject. As a result, our expressions allow for a direct comparison with the familiar 
Newtonian versions of Ohm's law, while identifying the relativistic corrections to them. We show, for example, that 
the relativistic analogue of the Biermann-battery effect has an additional contribution from the temporal pressure 
variations. Our main result is given in Eq. (j34p and applies to any multi-component fluid, relativistic or not, which 
means that it can be adapted to address a great variety of physical problems. With the general form of Ohm's in 
hand, our next step was to introduce a particular reference frame. Identifying our fundamental observers with the 
Eckart frame, allowed us to follow on the steps of the non-relativistic studies and therefore considerably simplify the 
mathematics. Then, by confining to two-fluid systems and assuming overall charge neutrality, we expressed Ohm's law 
in terms of the properties of the bulk. Finally, we closed our discussion by considered a number of applications. These 
included hot plasmas of two oppositely charged fluids, hot electron-positron mixtures, cold electron-proton systems 
and also the resistive and the ideal-MHD limits of our results. In each case, we have discussed the physics of the 
situation, identified the familiar effects, like the Biermann-battery and the Hall effects, and pointed to the relativistic 
corrections were applicable. 

The multi-fiuid description adopted in the present paper is essential in almost every study of small-scale astrophysical 
plasmas. The same approach is also necessary when looking into the nonlinear regime of galaxy formation, when the 
proto-structure has decoupled from the background expansion and collapses. Then, one can use our equations to 
investigate the evolution of proto-galactic magnetic fields, in particular their amplification and dissipation, within 
and also outside the MHD limit. Moreover, the use of the irreducible variables, assigns an unambiguous physical 
interpretation to every variable in our equations and helps to isolate the physical effects under consideration. For 
example, vorticity terms are related to turbulence and dynamo-like mechanisms, while those involving the shear 
describe shape distortions and can play an important role during galaxy formation. In addition, going beyond the 
cold-plasma limit, makes our equations suitable for studies of relativistic plasmas, like those in hot interstellar clouds 
and in accretion discs around compact stars. Applications of this sort will be the subject of future work. 

Acknowledgements: We thank Esteban Calzetta, Roy Maartens and Loukas Vlahos for useful discussions and 
comments. A.K. also acknowledges support from Projects PROPP-UESC 00220.1300.489 and 00220.1300.609. 

APPENDIX A: THE LANDAU FRAME 

As stated in § IV Al an alternative frame choice is that of the Landau (or energy) frame. Both the particle and the 
energy frames are physically equivalent and choosing one against the other depends on the particulars of the problem 
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in hand. The 4-velocities associated with the two frames are given by (j35p and p6p respectively, but rewriting Ohm's 
law relative to the Landau frame goes beyond the scope of this paper. We note, however, that the energy frame 
seems a less efficient choice, because both the interpretation of the resulting equations and their comparison to the 
non-relativistic expressions are less straightforward. Nevertheless, by projecting along m^, we can show that the 
same non-relativistic limit will be attained. Indeed, definitions (1351) . (1361) combine to 

„ r 1 



^[-n(,)A(,)+AAf,)#] . (Al) 



When the velocities of the individual species are small, we can neglect terms proportional to the heat flux. This means 
that for cold plasmas the above reduces to u^u^ ~ ^ / J2 ["'(j)A(j)] = ~1- Then, following (O, the Eckart 

and the Landau frames (and their associated expressions of Ohm's law) effectively coincide. 
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